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Abstract 

We investigate the possibility to draw conclusions on the sign of the spin- 
dependent gluon distribution, AG(x,Q 2 ), from existing polarized DIS data. The 
spin-dependent parton distributions Au v , Ad v , Au, Ad, As, and AG are constructed 
in the framework of a phenomenological procedure taking into account some as- 
sumptions on signs of valence and sea parton distributions motivated by 't Hooft's 
mechanism of quark-quark interaction induced by instantons. The axial gluon 
anomaly and data on integral quark contributions to the proton spin, Au, Ad, and 
As, are also taken into account. Predictions for the x- and Q 2 -dependencies of 
the polarized proton and neutron structure functions, g\ and g™, are compared to 
experimental data. It is shown that the neutron structure function, g™, is espe- 
cially sensitive to the sign of AG(x,Q 2 ). The results of our analysis supports the 
conclusion that this sign should be positive. 



1 Introduction 



Parton distributions in the nucleon are of universal nature, hence their parametrizations 
obtained from deep inelastic lepton-nucleon scattering can be utilized for simulations 
of processes outside lepton nucleon scattering; the polarized parton distributions are 
especially useful to predict the behaviour of pp interactions with polarized proton beams 
to facilitate future research programs at the RHIC, HERA and LHC colliders. 

Recent results deep inelastic lepton-nucleon scattering experiments at SLAC |], Q] and 
CERN [|J on spin- dependent structure functions for proton and deuteron targets, g\ and 
gf, stimulate the interest in determining the spin-dependent gluon and quark distribu- 
tions in a polarized nucleon. Since a complete solution of this problem is beyond the 
scope of perturbative QCD and there are still no sufficiently precise non-perturbative 
calculations available, the usual procedure is to fit numerous parametrizations of both 
spin-independent and spin-dependent parton distributions to the data. Up to now there 
is no unique solution; the results depend in one or the other way on the method used. 

Polarized parton distributions can be extracted in an indirect manner from doubly po- 
larized deep inelastic lepton-proton and lepton-deuteron scattering; the measurable ob- 
servables are the asymmetries A v and A d . The structure function g\ can be extracted 
from A p according to 

sfl(xm = A , (xm . -^gL_ Y (1) 

where additional information on the unpolarized structure function Q and on the ratio 
of longitudinal to transverse photon cross section R(x, Q 2 ) = <jl/&t are required. The 
deuteron structure function is determined in a similar way taking into account appropriate 
nuclear corrections. 

Since there is no practical way at present to directly extract polarized parton distributions 
from experimental data it is important to develop flexible procedures to construct these 
distributions incorporating relevant features of the data as well as reasonable constraints 
derived from our present theoretical understanding of the nucleon. 

At present there is no strong argument favouring a positive or negative sign of the spin- 
dependent gluon distribution, AG(x,Q 2 ). Several sets of spin-dependent parton distri- 
butions were constructed utilizing rather different approaches || 0, ||, U mostly assuming 
a positive sign of AG. Different parameter choices leading to a different behaviour of 



AG at x 1 (G t~ G I, G |> G j, G |< G j) were studied in |K|. Both positive 
and negative values of the sign of AG over a wide kinematical range 10~ 3 < x < 1 were 
considered in |§. A detailed NLO QCD analysis of the proton and deuteron data on gi 
was performed in || concluding that the size of the gluon distribution drives the pertur- 
bative evolution and, due to the fact that the SMC and E143 data were taken at different 
values of Q 2 , the polarized gluon distribution turned out to be large and positive. 

The aim of the present paper is to separate experimental observables being sufficiently 
sensitive to allow a determination of the sign of AG. As we shall show later, the neutron 
structure function g™ seems to be one of those observables. 



To construct the spin-dependent parton distributions a phenomenological method pro- 
posed in || is used. This method incorporates some constraints on shape and sign of 
parton distributions, it utilizes results on the quark contributions to the nucleon spin 
obtained in other analyses, and the effect of the axial anomaly is taken into account. We 
study the x and Q 2 dependence of g™(x,Q 2 ) for two different scenarios: AG > and 
AG < 0. The calculated predictions are compared to experimental data; a \ 2 criterion 
is chosen to judge in which of the two scenarios theoretical curves are better describing 
the experimental data on g™(x, Q 2 ). Eventually, the choice for a positive sign of AG will 
turn out to be the more likely one, i.e. the polarized structure function of the neutron 
will be shown to be sensitive to the sign of AG. 



2 Method 



The spin- dependent proton structure function g\ is expressed through spin-dependent 
parton distributions in a simple way 

9l(x, Q 2 ) = \- {^u(x, Q 2 ) + ±Ad{x, Q 2 ) + ±A~s(x, Q 2 )}, (2) 

where Aqf = qf — qj , and the qf are the probability distributions to find a quark having 
positive (+) or negative (— ) helicity relatively to positive proton helicity. The neu- 
tron structure function g±(x, Q 2 ) can be written in a similar form using the replacement 
Au <-> Ad. The valence distributions Aq v , Aq v are then obtained from Aq v = Ag — 2Aq. 
Since in this paper we shall use the spin-dependent parton distributions constructed in 
H we briefly describe in the following sections the main features of the applied method. 

2.1 Shape of Parton Distributions 

For the general form of a spin-dependent parton distribution Aqf we use 

Aqj = sign(qj) ■ x af ■ (1 — Xf f ■ qf, qf = u v , d v , u, d, s, G. (3) 

Here qf is the spin-independent parton distribution, are free parameters which are 

to be found by comparison with experimental data. From the restriction 

\Aq f \ < qf (4) 

follows that both probability distributions qf , qj as well as their sum qf = qf + qj 
need to be positive; moreover (3f should not be negative. To avoid the latter constraint 
we introduce a renormalised parton distribution qf = (1 — x)"i ■ qf. This leads to the 
following general form of a spin-dependent parton distribution 

Aqf = sign(qf) ■ x af ■ qf. (5) 

We note that since all presently available procedures to construct both spin-independent 
and spin-dependent distributions do imply fitting procedures and have consequently no 



unique solution. Hence we believe that at present it is recommended to incorporate gen- 
eral restrictions on Aqy like the one above; this makes it easier to develop flexible schemes 
to construct the helicity parton distributions and qj, separately. 



2.2 Signs of Parton Distributions 

Up to now there exists neither a running experiment to directly measure the polarized 
gluon distribution nor does the variety of indirect analyses give a unique result. Hence 
there exist no strong arguments on the sign of AG. Our approach will be to allow for 
both signs of AG and compare the quality of our model-dependent predictions to the 
experimental data. 

We note that a direct access to AG will be possible in future experiments. Utilizing 
polarized protons at RHIC for the (approved) experiments STAR and PHENIX [O] and, 



possibly, for the suggested internal polarized target experiment HERA-N |T2|] at HERA 



the measurement of AG seems feasible in the range 0.1 < x g i uon < 0.35. Also new doubly 



polarized lepton-nucleon scattering experiments proposed at CERN |L3| and suggested 
at SLAC |TJ], |T5| might contribute very valuable information on AG. 

For valence quark distributions the situation is much better defined; we take the sign of 
Au v as positive and that of Ad v as negative, respectively. This choice is motivated by the 
fact that the dominant configuration in the proton wave function is u('])u(])d(l), here 
the arrow denotes the quark spin direction. The same choice is made in most analyses of 



experimental data on quark contributions to the proton spin fll], 0, [|, fl6|j-fT9|j 



We assume for signs of Am (Ad) to be positive (negative). This is motivated by 't Hooft's 
mechanism p0[ for the spin configuration w(T) M (T)^(l) which determines the dynamics of 



quark helicity flips. The incoming left helicity quark q^ = (1 + 75)5/2 scattered from zero 
modes in the instanton field leads to an outgoing right helicity quark q^ = (1 — / y 5 )q/2. 
Effective Lagrangians are constructed in [ZT|; in the particular case of Nf = 2 flavours it 



can be written as 

L = J dp- n(p)(^ 2 p 3 ) 2 {u R u L d R d L [l + 1(1 - ^OAM + ( R ~ ( 6 ) 

Here p is the size of instanton, n(p) is the instanton density, = z/4- (7^7^ — 7^ 7 M ), and 
A a are matrixes for SU C (3) group. Once the left helicity quark scatters off an instanton 
it becomes a right helicity one and a quqn pair is created; the helicity of the sea quarks 
being opposite to that of the initial quark. In other words, the spin flip of the valence 
quarks u + and d~ determines the sign of the corresponding sea quark distributions - 
negative for Ad and positive for Au. A negative sign of As is in agreement with the 
arguments mentioned above and is supported by the results of several analyses of g{ data 



2.3 Inclusion of Axial Anomaly 



It was shown in [22| that the flavour-singlet axial current 



K = Yl 9/7//Y5?/ (7) 

f=u,d,s 



diverges at the quark level due to the one-loop triangle anomaly 



d » A = < *L.N r triF^F^}, (8) 

7T 

where F^ = e^ Pl F^ , F^ = d^A u - d v A^ + [A^A U ], A^ = A* ■ X a , a s is the strong 
coupling constant, and Nf is the number of flavours. The anomaly induces a mixing 
between the gluon and the flavour singlet axial current of quarks. For this reason, the 
helicity carried by each flavour undergoes a renormalization 



Aq f (x, Q 2 ) = Aq f (x, Q 2 ) - ^-L ■ AG(x, Q 2 ). (9) 

It was suggested in p2] that the axial anomaly might play a key role and modify the 
naive quark model predictions; hence parton distributions will presumably become much 
more sensitive to the sign of the polarized gluon distribution. Consequently, the spin- 
dependent structure functions g\ and g™ would become more sensitive to AG, as well. 



2.4 Integral Parton Contributions to the Proton Spin 

A further input required to our analysis is the total contribution of each quark species to 
the proton spin. We utilize the results of a recent analysis p3| of the structure functions 



g\ and gf from SMC and SLAC data incorporating 3 rd order pQCD corrections to the 
Bjorken sum rule. The relative quark contributions to the proton spin were determined 
as Am = 0.83 ± 0.03, Ad = -0.43 ± 0.03, As = -0.10 ± 0.03 at a renormalization scale 
Ql = 10 (GeV/c) 2 . Using these values and the definition 

f 1 Aq f (x, Ql)dx = Af, f = u,d,s (10) 

J 

the free parameters in the parametrization of our spin-dependent parton distribu- 

tions Any, Ad v , Au, Ad, As, AG were determined in M. 



3 Results and Discussion 

In fig. 1 (a,b) and 2 (a,b) the x-dependence of g\ and g x is shown for different parametriza- 
tions of parton distributions constructed with AG > (a) and AG < (b). The dashed, 
solid and dotted lines correspond to the parameters otf,(3f taken from Table 1-3 and 4-6 
of Ref . , respectively. 

From fig. 1 (a,b) is seen that all theoretical curves for the proton structure function g\ 
are in reasonable agreement with experimental data 0, |TB|, |T9"f , i.e. there seems to be 
no apparent sensitivity to the sign of AG. In contrast, from fig. 2 (a,b), displaying 
experimental data and theoretical curves for the neutron structure function xg™, one can 
deduce a certain dependence of the theoretical curves on the sign of AG in the range 
0.1 < x < 0.3. Hence there is some hope that xg™ exhibits a certain sensitivity to the 
sign of AG. 



Fig. 3 (a,b) shows the x- dependence of the proton structure function xg±(x, Q 2 ) at dif- 
ferent values of four-momentum transfer, Q 2 = 1,10,100 (GeV/c) 2 . The Q 2 behaviour 
of xg\ appears qualitatively different for AG > and AG < 0, respectively. In the first 
case the maximum of the curve is moved to lower x with increasing Q 2 , in the second one 
the position of the maximum is not affected. If AG > the prediction increases with Q 2 
for x < 0.01. If AG < 0, the prediction decreases with Q 2 over the full x-range. 

Fig. 4 (a,b) displays the ^-dependence of the neutron structure function xg™(x,Q 2 ) in 
the same fashion, i.e. for Q 2 = 1, 10 (GeV/c) 2 . If AG > the differences for different Q 2 
appear mainly at very low x-values and, in addition, at moderate x ~ 0.3. This sensitivity 
to the sign of AG is to weak for present experimental errors, however, it might be used 
later when more precise data on g™(x, Q 2 ) should become available. For AG < one 
observes a rather strong Q 2 -dependence at lower x-values and a somewhat characteristic 
dip at higher x, its position being almost independent on Q 2 . 

To be closer to the presently available Q 2 -values we show in fig. 5 the x-dependence of 
xg™(x, Q 2 ) at Q 2 = 1, 3, 5, 10 (GeV/c) 2 together with the presently available experimental 
data. (Due to the experimental errors the different ordinate is choosed in fig. 5(d) than 
in fig. 5(a,b,c).) The behaviour of xg™ on Q 2 is qualitatively and quantitatively different 
for the two scenarios AG > and AG < 0, especially at low Q 2 . Apparently, in the 
range x < 0.1 the experimental data on at Q 2 < 10 (GeV/c) 2 should be able to 
discriminate between positive and negative sign of the polarized gluon distribution. 

We apply a x 2 criterion to quantitatively distinguish between the two scenarios by com- 
paring our constructed parton distributions to the experimental data sets from SLAC and 
CERN |2|, |^, The obtained results are summarized in Table 1. There the references 
for experimental data on g™, the average Q 2 values, and the number of experimental 
points are shown in column 1, 2, and 3, respectively. The 'all' in col. 2 takes into account 
that each individual experimental point was measured at another average Q 2 , i.e. here 
the \ 2 is calculated using in the theoretical calculation the correct average Q 2 -value at 
each x-point. The corresponding kinematically accessible ranges are 1.1 5.2 (GeV/c) 2 
for E142 and 1.3 -j- 48.7 (GeV/c) 2 for SMC. This method seems to us the closest pos- 
sible description of the data by a theoretical calculation, hence we expect the x 2 values 
obtained for the 'all' comparison to yield the best possible separation. 



Experiment 


<Q 2 > 
(GeV/c) 2 


data 
points 


X 2 1 ndf 
AG > 


X 2 1 ndf 
AG < 


E142 [18| 


2 


8 


1.20 


2.05 


E143 [| 




3 


9 


0.89 


1.41 


SMC § 


10 


12 


1.28 


1.63 


HERMES |24 


3 


8 


0.86 


1.20 


E142 [18| 


all 


8 


1.45 


2.30 


SMC | 




all 


12 


1.35 


2.41 



Table 1. x 2 comparison between theoretical predictions, calculated for the two scenarios 
AG > and AG < 0, and experimental data on (x, Q 2 ). 



From table 1 one can see that for every data set the \ 2 P er degree of freedom is signif- 
icantly better in the case AG > compared to the case AG < 0. These results can be 
considered as clear quantitative evidence that the case AG > is the more likely scenario 
compared to the case AG < 0. 

We note that our result supports the conclusion on a positive sign of AG > obtained 
recently by a NLO QCD fit to g\ proton and deuteron data ||. 

Finally we present in table 2 our results for the integral quark - AS - and gluon - Ag - 
contributions to the proton spin calculated with the low limit x m i n = 10~ 3 . Whereas in 
the more likely scenario AG > the quark contribution AS appears to be almost stable 
with Q 2 it drops by almost a factor of 2 when increasing Q 2 from 3 to 10 (GeV/c) 2 in 
the less likely case AG < 0. In both scenarios Ag rises by about 10% within the same 
Q 2 range. 



Qo 

(GeV/c) 2 


AS 


Ag 


AG > 


AG < 


AG > 


AG < 


3 


0.290 


0.520 


1.78 


-3.01 


5 


0.293 


0.420 


1.86 


-3.20 


10 


0.298 


0.296 


1.95 


-3.41 



Table 2. Integral quark - AS - and gluon - Ag - contributions to the proton spin, cal- 
culated from the constructed polarized parton distribution functions for the two scenarios 
AG > and AG < 0. 



4 Conclusions 

The possibility to draw conclusions on a positive or negative sign of the polarized gluon 
distribution AG (a;, Q 2 ) was studied using a phenomenological procedure to construct 
spin-dependent parton distributions. The method includes some constraints on the signs 
of valence and sea quark distributions, takes into account the axial gluon anomaly and 
utilizes results on integral contributions to the nucleon spin, Am, Ad, As. Investigating 
the x- and Q 2 -dependencies of the structure functions g\ and g™ constructed by this 
method we introduce a x 2 criterion to discriminate between the two scenarios obtained 
for AG > and AG < 0, respectively. The neutron structure function turned out to be 
sufficiently sensitive to the sign of AG(x, Q 2 ), even at the present level of only moderate 
experimental errors. The results of our analysis strongly support the conclusion that the 
sign of AG(x, Q 2 ) is positive. New data on the neutron structure function g™ from the 
latest SLAC experiments and from HERMES at DESY will undoubtedly allow to draw 
a more definite conclusion on the sign of the polarized gluon distribution. 
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Figure 1. Deep-inelastic proton structure function xg\(x,Q 2 ). Experimental data: * 
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a) - AG > and (b) - AG < at Q 2 = 
Parametrizations of parton distributions: - - , - - , - - are taken 
from Tables 1-3 and Tables 4-6 M, respectively. 
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Figure 2. Deep-inelastic neutron structure function xg™(x, Q ). Experimental data: • 
- 0. Theoretical curves: (a) - AG > and (b) - AG < at Q 2 = 10 (GeV/c) 2 . 



Parametrizations of parton distributions: - - , - - , - - are taken from Tables 1-3 
and Tables 4-6 M, respectively. 
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Figure 3. Deep-inelastic proton structure function xgi(x,Q 2 ). Experimental data: * 
- 0, • - 0, o - HJ. Theoretical curves: (a) - AG > 0, (b) - AG < and 
1 (GeV/c) 2 , - - 10 (GeV/c) 2 , - - 100 (Ge1//c) 2 . Parametrizations of parton 

distributions are taken from Tables 2 and 5 ||. 
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Figure 4. Deep-inelastic neutron structure function xg™. Experimental data: 
Theoretical curves: (a) - AG > 0, (b) - AG < and - - 1 {GeV/cf 
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(GeV/c) 2 . Parametrizations of parton distributions are taken from Tables 2 and 5 
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Figure 5. Deep-inelastic neutron structure function xg™(x,Q 2 
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Experimental data: * 
H|. Theoretical curves: - - - AG > 0, — - AG < at 



Q 2 = 1, 3, 5, 10 (GeV/c) 2 . Parametrizations of parton distributions are taken from Tables 
2 and 5 [|. 



0.03 




c) 




d) 



Figure 5. Continued 



